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I. INTRODUCTION 

Extended one-dimensional (ID) models are excellent 
approximations for the electronic structure of some crys- 
tals, either inorganic or organic. Quite separately, ID 
models have interesting theoretical and thermodynamic 
properties. In addition to exact results, approximate 
methods have been widely applied to and tested on ID 
models. Two major recent developments are the density 
matrix renormalization group (DMRG) and field theory. 
The two methods are complementary in principle, and 
both have been applied to the zigzag spin-1/2 chain that 
is the subject of this paper. In practice, however, field 
theory deals with small energy gaps or long correlations 
lengths near quantum phase transitions that may be 
beyond the accuracy of numerical methods, a point 
often made for Kosterlitz-Thouless transitions. The two 
approaches to extended ID systems are quite different. 
DMRG is a versatile numerical technique for growing an 
extended chain from a finite one. It provides a complete 
approximate description of the ground state (gs) or 
other properties. When multiple DMRG schemes are 
possible, the most accurate one is readily identified. 
Field theory is an analytical approach based on a 
continuum approximation, or an effective Hamiltonian, 
to a discrete ID model. It targets critical phenomena at 
quantum phase transitions. A ID model may support 
multiple field theories among which it may be difficult 
to choose. 

In this paper, we present a modified DMRG algorithm 



to the zigzag spin-1/2 chain with frustrated antiferro- 
magnetic (AF) exchange Ji > and J2 > between first 
and second neighbors. The Hamiltonian of this familiar 
ID spin system is 

H{x) = [(1 - x)Sn ■ Sn+1 + xSn ■ Sn+2\- (1) 

n 

We consider the interval < x < 1 and set the total 
exchange J = 1 as the unit of energy. The x = limit 
is a linear Heisenberg antiferromagnet (HAF) with many 
known exact properties fl*! and many physical realiza- 
tions. The X K, 1 limit corresponds to two HAFs, one on 
each sublattice, and is the zigzag chain sketched in Fig. 
[T] Small Ji for x < 1 or x > 1 describes an interchain 
exchange that is AF or ferromagnetic (F), respectively, 
and is frustrated because each spin is equally coupled 
to two neighbors of the other sublattice. The modified 
algorithm improves the accuracy for x > 0.5 (J2/J1 > 1). 

The spin chain H[x) has been extensively studied, es- 
pecially in the x = limit [l| that Bethe [2l| and Hulthen 
3] solved long ago. Majumdar and Ghosh (MG) found 
J] a simple exact gs at xmg = 1/3 (J2 = >/i/2). The 
gs is a doubly degenerate bond order wave (BOW) with 
broken inversion symmetry at sites. The fluid-dimer 
transition with increasing J2/J1 marks the opening of 
a magnetic gap A(x) between the singlet gs and the 
lowest triplet state. The first field theoretic treatment 
j5] placed the critical ratio of J2/J1 at 1/3; subsequent 
analysis returned ^] J2IJ1 = 1/6 and finally 7] « 1/4. 
Okamoto and Nomura Q obtained the accepted value, 
J2/J1 — 0.2411 or Xc = 0.1943 in our notation, using 
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FIG. 1: Zigzag spm-1/2 chain with AF exchange Ji and J2 
between first and second neighbors, respectively. Ji = gives 
two Heisenberg chains of 2n spins with exchange J2 



exact results up to iV = 24 sites, extrapolation and field 
theory. 

Classical spins in H{x) lead to spiral phases for large 
J2I J\ when adjacent spins are nearly orthogonal. The gs 
energy per site for classical spins with angle Q between 
neighbors is 



(2) 



Minimization with respect to G leads to cosB = — J1/4J2 
and Q = tt/2 +_y for large J2/J1. The spiral phases of 
quantum spins |9l4ll| are another area of interest, as are 
the structure factor S{q) as a function of J2/J1 and 
the crossover [l^ from a singlet to a ferromagnetic gs at 
Ji = — 4J2. There are possible physical realizations [13] 
of H{x), most with AF exchanges Ji, J2 and a few with 
F exchange Ji. 

White and Affleck (WA) studied [U] the BOW phase 
with J2 / Ji beyond the MG point by a combination of 
DMRG and field theory. Numerical issues limited DMRG 
to J2/J1 — 2.0 for A and to 2.5 for the order parame- 
ter. The modified algorithm is accurate up to J2/J1 = 4. 
WA concluded that the BOW phase extends to a; = 1 
( Ji = 0). Itoi and Qin (IQ) presented [l^ a more elabo- 
rate field theory for large J2/J1. The present work was 
motivated in part by the contrasting results of IQ and 
WA. According to IQ, the spin correlation length diverges 
as [H 

e(Ji, J2) «ea;p(c(|Ji|/J2)~'/'), (3) 

where c is a constant. The WA expression for ^ has 
exponent —1 instead of —2/3 and is limited to Ji > 0. 
The order parameter of the BOW phase is 



B{x) = {Sn ■ Sn+l) - {Sn ■ Sn-l)- 



(4) 



B{x) is the gs amplitude of the BOW for x > Xc- WA 
call it "dimerization" , a term that we reserve (l6i| for 
structurally dimerized systems such as polyacetylene or 
ion-radical salts or spin chains. Broken inversion symme- 
try in a BOW phase is electronic dimerization in a regular 
array. Both WA and IQ support their ^{x) with the same 
(limited) DMRG resuhs [14] for B{x) and A (a;). 

Since B{x) and A(a;) are proportional to l/^(x), 
DMRG for the BOW amplitude or the magnetic gap can 



be compared to field theory as 

lnB{x) « lnA{x) « -c( Ji/ J2)~^/^ (5) 

with J2/J1 = x/{l — x). The numerical problem is to 
evaluate exponentially small quantities at large J2/Ji- 
The two computations are independent, since A(a;) re- 
quires the triplet state while B{x) does not. DMRG di- 
rectly yields approximate spin correlations functions in 
the gs 



C{p) = {Sn ■ Sn+p), 



(6) 



and the wavelength \{x) of a spiral phase, if present. As 
noted by WA [IJ], B{x) and A(a;) are well-defined quan- 
tities whereas £,{x) requires an unknown fitting function 
in addition to C{p). The order parameter of the spiral 
phase is the twist angle x{x) below Eq. ^ that is related 
Ini El to the BOW phase as 



e{x) 



= X{x) = 



4^ 



(7) 



x(a;) — 2t:/X{x) has been approximated by a coupled- 
cluster expansion J9| and by twisted boundary conditions 
in finite systems [llj . 

A spiral phase of H{x) has been analyzed [l^l in the 
classical limit of an infinite spin at each site in terms of a 
nonlinear cr-model that involves a 3 x 3 orthogonal matrix. 
That field theory does not produce a BOW, however, and 
is not powerful enough to yield scaling results for X{x) or 
A(a;). On the other hand, field theories [13: [iB] based on 
bosonization do not predict a parameter range in which 
a spiral phase should appear. Indeed, there is no com- 
pelling field theoretic reason that necessarily relates the 
spiral and BOW phase. They have different order pa- 
rameters and different symmetries, a discrete symmetry 
for translation by one site in the BOW phase and a con- 
tinuous rotational symmetry for the spiral phase. The 
BOW extends from Q = 0.1943 to [3 a; = 1, while 
the range of a spiral phase is [13] from xmg — 1/3 to 
a; = 1. The richness of the zigzag chain at large J2/J1 
makes it ideal for a critical discussion of DMRG accuracy 
and comparisons to field theory. 

Section II describes the modified DMRG algorithm in 
which four rather than two spins are added per step. The 
accuracy improves modestly at a; = and dramatically 
for a; > 2/3 where the second-neighbor J2 dominates. 
Adding four spins when Ji is small amounts to increas- 
ing two weakly-coupled chains by two spins each, just as 
adding two spins does at a; = in conventional DMRG. 
We present results in Section III for B{x) and A (a;) up 
to X = 0.8 (J2/J1 = 4) and for xi^) ^-^d A(x) up to 
X = 0.75 (J2/J1 — 3), the practical limit in chains of 
N < 1000 spins with open boundary conditions (OBC). 
Our results agree with the IQ expression in Eq. ^ with 
c — 2.90 ± 0.10 for all four quantities. We also com- 
pute the structure factor S{q) and its maximum q* that 
yields an independent estimate of the twist angle x(a;). 
We comment in Section IV on the status of comparisons 
between DMRG and field theory for the zigzag chain at 
large J2 / Ji . 



II. MODIFIED DMRG ALGORITHM 

DMRG is among the most accurate numerical tech- 
niques for solving extended ID quantum cell models [isl - 
[2H . Conventional DMRG algorithms start with four sites 
and grow an extended chain by adding two sites in the 
middle, treating the left and right half-blocks as system 
and environment by turns [3, [I^ . The accuracy of this 
method decreases for long-range (beyond first neighbor) 
interactions since we encounter bonds between old sites 
of the same system block. Long-range interactions in 
conventional DMRG couple sites at every step whose op- 
erators have undergone an unequal number of renormal- 
izations. The spin chain H{x) in Eq. ([1]) has second- 
neighbor J2 between sites introduced on successive steps. 
The decrease in accuracy becomes significant when J2 is 
large because site operators involved in J2 are renormal- 
ized twice while the Ji operators are renormalized only 
once. A remedy is to add sites on every step that encom- 
pass the full range of interactions. 

Accordingly, we modified the DMRG algorithm for 
H{x) to add two new sites per half block instead of one, 
as shown schematically in Fig. [2] The system starts with 
4 spins and grows to = 4ri in iV — 1 steps. Large J2/J1 
leads to weakly coupled chains of 2n sites in Fig. (TJ each 
with a singlet gs when Ji — 0. The Fock space dimen- 
sionality of each block increases as 2(2s -I- l)m, or as 4m 
for s = 1/2 sites, which is comparable to fermionic sys- 
tems. We find that keeping m = 150 eigenvectors of the 
density matrix is sufficient for good accuracy. The trun- 
cation error in the sum of the eigenvalues of the density 
matrix is less than 10~^ in the worst case, and increasing 
m changes the energy only in 5*^ or 6"* decimal place, in 
units of J. For more accurate spin correlation functions 
C{p) and order parameter B{x), we used finite DMRG 
calculations on every fourth steps [l^. B{x) is calculated 
using the middle bonds of the chain and is accurate up 
to 5-6 decimal place, but is subject to finite-size effects 
of order 1/A^ discussed below. 

We compare the modified algorithm with four sites 
added per step to conventional DMRG for the gs energy 
and the gap magnetic A. We use the infinite DMRG al- 
gorithm with m = 200 in each case. Since DMRG targets 
the lowest state in each Ms sector, the lowest permissi- 
ble total spin state has the best energy in an AF model. 
Conventional DMRG is most accurate aX x = 0(J2 = 0) 
where there is only nearest-neighbor exchange. Neverthe- 
less, as shown by the inset in Fig. [3J the new method im- 
proves the gs energy slightly and the triplet energy con- 
siderably. Note that the inset energy scale is 100 times 
finer that of the main figure. We attribute better per- 
formance to (i) the absence of old-old bonds within the 
same block in the new scheme, and (ii) increased num- 
ber of new-new bonds (3 at a; = 0) compared to new-old 
bonds (2 at a; = 0) when four new sites are added at 
each step. The conventional ratio is 1:2 at x = 0. The 
accuracy of the new method at x = is about 10~^ for 
the singlet and 10~^ for the triplet. It runs smoothly for 
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FIG. 2: DMRG scheme with four new sites added per step. 
Primed and unprimed indices are sites of the left and right 
blocks, respectively. Open circles represent new sites and 
closed circles, old sites. Solid lines represent J2, dashed lines 
Ji. 
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FIG. 3: Energy difference per site between the new {E4) and 
conventional (£'2) DMRG for the singlet gs and the lowest 
triplet at X = 0.70 and x = (inset) for chains of A'^ sites. 



X > 2/3, in contrast to numerical difficulties [T4| of con- 
ventional DMRG at X > 1/2. The estimated accuracy for 
X > 0.5 is 10^^ for the gs and lO^'^ for the triplet. We 
also studied chains with J2 > and F exchange Ji < 
in terms of Ji, J2 rather than a; > 1 in Eq. ([1]). 

Figure S] shows the size dependence of B(0.8) and 
A(0.8) at J2/J1 = 4.0. These are the smallest B and 
A that are accurate with the present DMRG. We varied 
TO to look for jumps in B{x), such as those in Fig. [S] 
of ref. 14 at J2/J1 = 2.5, but found only the smooth 
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FIG. 4: BOW order parameter B(0.80) vs 1/N for finite and FIG. 5: Order parameter B and magnetic gap A as a function 

infinite DMRG and (inset) magnetic gap A(0.80) vs l/N. of (Ji/J2)"^/^. The fitted lines are l/£,{x) in Eq. (O with 

Infinite DMRG is carried out with m = 300 for N > 200 and c = 2.94 for A and 2.86 for B. 
m = f 50 for iV < 200; finite DMRG has m = 200. 



behavior shown. In the four-site algorithm, B(x) is hn- 
ear in 1/7V for large N. Finite DMRG procedure with 
m = 200 and N between 100 and 200 sites returns 
5(0.8) = 0.0071, as shown in Fig. H The infinite al- 
gorithm with variable m, and 200 < N < 430 leads to 
extrapolated 5(0.8) 0.0066. The inset of Fig. Ushows 
the 1/N dependence of A(0.8) using finite DMRG with 
four spins added per step. The extrapolated gap is 0.002. 
Similar extrapolation at a; = give A w 0.001, close to 
the exact A = 0. 



III. THE J2/J1 > 1 REGIME OF H{x) 

Larger J2/J1 is accessible with the improved DMRG 
algorithm. All results below are for m — 200 and OBC 
for N — 800 sites, as discussed in Section II. The order 
parameter B{x) in Eq. ^ and the magnetic gap A(a;) 
from the gs to the lowest triplet provide direct compar- 
ison to field theory of the BOW phase. Both B{x) and 
A(a;) go as l/^(x), where ^{x) is the correlation length 
in Eq. ([3]). As seen in Fig. [SJ the IQ exponent of —2/3 
fits the DMRG results remarkably well up to J2/J1 =4 
{x = 0.8) with c — 2.90. The B{x) fit covers almost two 
orders of magnitude and extends to J2 = Ji- 

The J2/J1 = 1 point for A deviates upward from the 
line in Fig. [S] IQ used DMRG [13 for A(x) in the 
interval 0.6 < J2/ Ji < 2 to support ^ in Eq. ^ with c = 
3.66. This is not correct, and the one-loop approximation 
does not extend down to J2/J1 ~ 1. WA used DMRG 
for B(x) up to J2/J1 2.5 to support ^ with exponent 
— 1 instead of —2/3 in Eq. Their expression fails 

at larger J2/J1. DMRG with 0.6 < J2/J1 < 2.5 is not 



appropriate for the BOW phase at large J2/ Ji- Although 
Fig. O covers more than an order of magnitude in A and 
almost two for B, there is no assurance that J2I J\ = 4 
is large enough. On the other hand, if J2/J1 ~ 4 is 
not "large" , numerical comparison with field theory will 
indeed be difficult. 

The accuracy of B{x) is limited by finite-size effects 
for OBC and N — 800 sites. We illustrate with an un- 
correlated example. A half-filled Hiickel or tight-binding 
chain of N sites has bond orders [l^l 

N/2 

Pm — '^'^^Ck,mCk,m+l, (8) 
fc=l 

with m = 1,2,..., — 1. The coefficient Ck,m at site m 
of the filled orbital k is 

/ 2 ivkm 
Ch m = \ sin . (9) 

' V Af + 1 A^+ 1 

The geometrical series for pm is summed for finite N . The 
difference between Pm/2 of the central bond and that of 
either neig hbor is -2(-1)^/VA^ for large A^. The bond 
order Pm/2 is less than the band limit of 2/ti for N — An 
and greater than 2/7r for A^ = An + 2, just as expected 
for partial single and double bonds at the center of linear 
polyenes with evenly spaced C atoms. Since OBC break 
inversion symmetry at sites, this elementary example has 
implications for any OBC simulation of BOW systems. 
In any case, the exponential decrease of B with J2/J1 is 
soon overwhelmed by 1 /N corrections that limit DMRG 
with N « 1000. Finite-size corrections to A or other low- 
energy excitations also go as « 1/N and place similar 
limits on the accuracy of exponentially small gaps. 
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FIG. 6: Spin correlation functions C(p) on the same sublattice 
(open symbols, even p) and on opposite sublattices (closed 
symbols, odd p) for A'^ = 800 sites. The scaling of C{p) and 
choice of are discussed in text. 



A DMRG calculation returns all gs spin correlations 
functions C{p) in Eq. OBC implies that C{n,p) 

depends on the site index n as well as the separation p. 
It is customary to take sites n and n + p in the central 
part of the chain. C{p) between sites on one sublattice 
in Fig. [T] has even p, while C (p) between sublattices has 
odd p. Fig. [S] shows C{p) in spiral phases at a; = 0.65 
in the top panel and x = 0.675 in the bottom panel. 
The wavelength A(x) of the spiral phase appears directly 
provided that there are two nodes to specify A/2. DMRG 
with N = 800 sites yields A only up to x = 0.75. The 
scale factor p^/^ea;p(p/^) follows WA, who considered 
even p and chose ^ to make the amplitudes in Fig. [6] as 
equal as possible. The same ^ holds for odd p. This 
procedure minimally requires two maxima and hence is 
also limited to x — 0.75(J2/Ji — 3). 

We obtained A(x) and ^(a;) from C(p) results with N = 
800 sites. DMRG is unbiased in the sense that neither a 
spiral nor a BOW phase is assumed. Fig. [7] shows how 
ln^{x) and ln\{x) increase with Jij J\. The IQ exponent 
of —2/3 in Eq. ([3]) fits reasonably well over a smaller 
range of J2/J1 with c = 3.03 for A. The ^ exponent 
is consistent with the more accurate c = 2.90 in Fig. 
[5] for B{x) and A(a;) over a wider range. The scaling 
form of ^(x) in the BOW phase and \{x) in the spiral 
phase are almost identical. According to Eq. ([7]), the 
product \{x)B{x) or of \{x)^{x) should be constant, 
independent of x for large J2 / Ji . The calculated points 
in Figs. [S] and [7] between J2/J1 — 1.3 and 3.0 yield 
\{x)B{x) PS 9 ± 3 and A(a;)A(a;) w 5.6 ± 2. Since neither 
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FIG. 7: Spin correlation length ^ of the BOW phase and 
the wavelength \(x) of the spiral phase as functions of 
{J\l J^Y'^^'^ , based on the C(p) results on Fig. [6] 



is monotonic in J1I J\i our results are weakly consistent 
with constant \[x)l^{x). Higher accuracy is needed to 
test Eq. dZl). 

The spiral phase of H[x) has been modeled [9l-[ll| in 
terms of the twist angle x > for AF exchange that is 
defined below Eq. ^ for classical spins. The inverse 
relation between \{x) and ^{x) in Eq. (O has been pro- 
posed [niQ for large A (a;) or small x when the discrete 
nature of the spin chain is irrelevant. It follows that 

C(2p) oc cos{2pd) ^ {-\fcos{2px), 
C(2p-|-1) oc cos[(2p+ 1)61] = -(-l)Pcos(2px)-(10) 

Even and odd C{p) are not quite out of phase in Fig. [51 
in agreement with Eq. ([TU]). The nodes of C(2r) occur 
at 2rx = (n + l/2)7r, while those of C(2r + 1) are at 
(2r+l)x — mr. The angle x(a:) decreases with increasing 
J2/J1 and has been studied by other techniques [9, 11]. 
Independent evaluation of x(a;) — 2t:/X{x) provides a 
consistency check for direct DMRG results for A(a;) in 
Fig. [6l Such consistency is different from the common 
scaling of S,{x) and A(x) discussed above. 

AligiaetaL [ll| obtained x(a;) using twisted boundary 
conditions in Eq. ((TJ and exact results up to = 24. 
Bursill et al. Q presented several approximation schemes 
for x{^): one of which is based on the peak q* of the 
structure factor S{q). The spin-1/2 structure factor for 
a system with periodic boundary conditions is 

^(«) = ^ E cip)^M^qp) = l + J2 ^cip)cos{qp)in) 

np p—1 

where C{p) are spin correlation functions in Eq. ([SJ. 
Inversion symmetry is restored in a BOW phase by taking 
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FIG. 8: Derivative S'{q) of the structure factor. The inset 
shows S{q*) = for a; = 0.625, 0.65 and 0.675. The twist 
angle x is 9* — 7!"/2 in radians. The x = 0.700 result is un- 
physical. 



a linear combination of the degenerate gs. The MG point 
at a; = 1/3 has short-range correlations, known exactly, 
leading to Smg{q) — 3{l — cosq)/4: and a broad maximum 
at q* — TT. The maximmn value S{q*) is obtained using 
the derivative 



dSjq) 
dq 



^2pC{p)sin{qp). 



(12) 



Eq. ([T2|) shows that q* is sensitive to long-range spin 
correlation functions. We again use C{p) from the central 
part of the chain. C{p) refers to n = N/2 = 400 in Eq. 
(|6]) and the sum is fromp = 1 to iV/2 — 10, or 10 sites from 
chain end. The resulting S'{q) are shown in Fig. [S] The 
inset magnifies the S'{q*) = region for the indicated 
values of x. As ^(x) increases and correlations become 
long ranged, large p must be retained in the sum and the 
inherent 1/iV limitations of OBC are again encountered. 
Although q* — tt/2 + xi^) — J" 7r/2 with increasing x as 
expected, the condition S'{q*) = has limited value in 
the crucial region of small x- The point q* — tt/2 occurs 
at Ji = that separates the AF regime with Ji > 
and q* > tt/2 from the F regime with Ji < and q* < 
tt/2. We underestimate q* for x = 0.70, which is clearly 
unphysical, and hence overestimate A based on S{q), but 
the twist angle and wavelength are consistent for x < 
0.65. 

Aligia et al. llj emphasize that twisted boundary con- 
ditions extend xi^) to much larger J2/ Ji ~ 30. They re- 
port reasonable agreement with WA IJ] and with Bursill 
et al. up to J2/ Ji = 2.5, where our results are simi- 
lar. But A(0.8) ~ 1300 estimated from J2/J1 = 4 in their 
Fig. 4 is about 7 times smaller than the extrapolation of 



A (a;) in Fig. [T) Moreover, their [iTj] asymptotic regime 
starts at J2/J1 — 15 where their ^(x) has the WA form 
with exponent —1 in Eq. ([3]). The stronger presumed de- 
crease of l/£^{x) in the spiral phase as x — ?> 1 would lose 
out to the weaker singularity of the BOW phase. Twisted 
boundary conditions up to A'^ w 24 do not to give reliable 



11| A(x), however, and no B{x) results were presented. 



IV. DISCUSSION 

We obtained more accurate results for the frustrated 
spin chain H{x) in Eq. ([1]) with J2/J1 > 1 by modifying 
the DMRG algorithm to add four sites per step instead 
of two. The order parameter B{x) in Eq. ^ is limited 
by 1 /N corrections in systems with open boundary con- 
ditions (OBC). The accuracy of the magnetic gap A(x) 
to the lowest triplet is estimated by comparison to exact 
results in the fluid phase with x < Xc = 0.1943. As seen 
in Fig. [5l we find an exponential decrease of B{x) and 
A(x) up to J2/J1 — 4 that follows the IQ [lB| correla- 
tion function in Eq. (jS]) for almost two decades. DMRG 
automatically yields the spin correlation functions C{p) 
in Eq. (jG)) and a spiral phase with wavelength \{x) in 
Fig. ini Following WA [1J|, the correlation length ^(x) is 
extracted from amplitudes in the spiral phase. Exponen- 
tially increasing A(x) and ^(x) in Fig. [7]again follows the 
IQ expression in Eq. ([3]) , albeit over a narrower range up 
to X = 0.75(J2/Ji = 3) set by numerical considerations. 
The maximum q* of the structure factor in Fig. [5] is an 
independent estimate of twist angle xi^) = 2tt/\{x) of 
the spiral phase in Eq. ([2]). We find that q* has limited 
accuracy for our C{p) for x > 0.65. 

Detailed comparison with theory is made possible by 
multiple studies of the BOW [11 lH and spiral [1, [El 
phases of H{x) with J2/J1 > 1- More generally, we won- 
dered whether DMRG is capable of confirming the small 
gaps or long correlation lengths predicted by field theory. 
Our results to J2/J1 = 4 clearly favor the IQ expression 
[Tsl for i (x) in Eq. ([3]) while just as clearly ruling out 
their fit [15l| for A(x). Greater accuracy is needed for 
meaningful comparisons. The modified algorithm yields 
multiple and reasonably consistent comparisons up to 
J2/J1 = 4. 

The modified algorithm runs smoothly for J2 > and 
Ji < 0. The IQ expression for ^(x) in Eq. ^ does not 
depend on the sign of Ji . Wc find finite gaps A on the F 
side that, however, are less than our estimated numerical 
accuracy. Still higher accuracy is needed for A on the F 
side. We can definitely say, however, that the constant 
c ~ 2.9 for ^(x) on the AF side is different from that on 
the F side. In view of small A, Itoi and Qin discuss [lH| 
the spin wave velocity of the singlet or triplet and present 
conventional DMRG results for A^A vs. l/N in Figs. [3] 
and ID of Ref. 0. The appearance of a nonsinglet gs at 
Ji — 2J2 contradicts the exact result of Dmitriev et 
al. [l3], that the singlet/ferromagnetic phase boundary 
of the zigzag chain is at Ji = — 4J2- Field theory on the 
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F side is numerically untested so far. 

DMRG accounts naturally for coexisting BOW and 
spiral phases with onsets at Xc = 0.1943 and xmg = 
1/3, respectively, but cannot say where they terminate. 
Bosonization field theories [l3, ISl have a BOW phase 
but not a spiral phase, while field theory [lo| or other ap- 
proaches [9, TT] to the spiral phase do not yield a BOW. 
Eq. ([7]) is an assumed 11, 14] relation between the twist 
angle x{^) of the spiral phase and the correlation length 
^(x) of the BOW phase. The scaling of \{x) and ^(a;) 
in Fig. [7]is almost the same, and the products X{x)B{x) 
and A(x)A(a;) are roughly constant, but greater accuracy 
is needed to confirm that A(a:;) and S^(x) are indeed pro- 
portional. It may be interesting in the future to study 
whether similar scaling is special to spin- 1/2 or holds also 
for higher spin. 

Field theory is a continuum approximation. Since 
solid-state models are discrete, field theory becomes ac- 
curate when ^{x) exceeds 5-10 lattice constants. This is 
well documented for solitons in the SSH model [23] and 
its continuum version [23|. As shown in Fig. [71 ^(x) > 10 



requires Jij J\ > 1 and our DMRG extends to ^ w 300. 
We do not consider the discreteness of the lattice to be 
important. 

It is a well-recognized numerical challenge, to obtain 
exponentially small energy gaps or exponentially long 
correlation lengths near a quantum phase transition. 
Impressive gains in numerical accuracy are required to 
be modestly closer to the critical point. The modified 
DMRG algorithm for H{x) extends accurate results to 
J2/J1 — 4 and clearly favors the correlation function 
^{x) in Eq. ^ proposed by Itoi and Qin [l^. There 
are open questions such as whether J2/J1 = 4 is in the 
asymptotic limit or the relation between BOW and spi- 
ral phases. Convincing comparison between field theory 
and numerical methods are in fact demanding as we have 
illustrated for the zigzag spin-1/2 chain. 

Acknowledgments. We gratefully acknowledge par- 
tial support for work at Princeton by the National Sci- 
ence Foundation under the MRSEC program (DMR- 
0819860). SR thanks DST India for funding through 
SR/Sl/IC-08/2008 and JC Bose fellowship. 



[1] D. C. Johnston, R. K. Kremer, M. Troyer, X. Wang, 
A. Klumpfer, S. L. Bud'ko, A. F. Panchula and P. C. 
Canfield, Phys. Rev. B 61, 9558 (2000), and references 
therein. 

[2] H. Bethe, Z. Phys. 71, 205 (1931). 

[3] L. Hulthen, Arkiv. Mat. Astron. Fysik. 26A, No. 11 
(1938). 

[4] C. K. Majumdar and D. K. Ghosh, J. Math. Phys. 10, 
1399 (1969). 

[5] F. D. M. Haldane, Phys. Rev. B 25, 4925 (1982). 
[6] K. Kuboki and H. Fukuyama, J. Phys. Soc. Japan 56, 
3126 (1987). 

[7] I. Affleck, D. Gepner, H. J. Schultz and T. Ziman, J. 

Phys. A 22, 511 (1989). 
[8] K. Okamoto and K. Nomura, Phys. Lett. A 169, 433 

(1992). 

[9] R. Bursill, G. A. Gehring, D. J. J. Farnell, J. B. Parkin- 
son, T. Xiang and C. Zeng, J. Phys: Condens. Matter 7, 
8605 (1995). 

[10] S. Rao and D. Sen, Nucl. Phys. B 424, 547 (1994); S. 

Allen and D. Senechal, Phys. Rev. B 51, 6394 (1995). 
[11] A. A. Aligia, C. D. Batista and F. H. L. Essler, Phys. 

Rev. B 62, 3259 (2000). 
[12] D. V. Dmitriev, V. Ya. Krivnov and A. A. Ovchinnikov, 



Phys. Rev. B 56, 5985 (1997). 
[13] M. Hase, H. Kuroe, K. Ozawa, O. Suzuki, H. Kitazawa, 

G. Kido and T. Sekine, Phys. Rev. B 70, 104426 (2004). 
[14] S. R. White and I. Affleck, Phys. Rev. B 54, 9862 (1996). 
[15] C. Itoi and S. Qin, Phys. Rev. B 63, 224423 (2001). 
[16] M. Kumar, Z. G. Soos and S. Ramasesha, Phys. Rev. B 

(in press). 

[17] T. Tonegawa and I. Harada, J. Phys. Soc. Japan 56, 2153 

(1987); R. Ghitra, S. Pati, H. R. Krishnamurthy, D. Sen 

and S. Ramasesha, Phys. Rev. B 52, 6581 (1995). 
[18] S. R. White, Phys. Rev. Lett. 69, 2863 (1992). 
[19] S. R. White, Phys. Rev. B 48, 10345 (1993) 
[20] U. Schollwock, Rev. Mod. Phys. 77, 259 (2005) 
[21] K. Hallberg, New Trends in Density Matrix Renormal- 

ization, Advances in Phys. 55, 477 (2006). 
[22] C. A. Goulson, Proc. R. Soc. London, A 169, 413 (1939); 

Z. G. Soos, S. Ramasesha, D. S. Galvao and S. Etemad, 

Phys. Rev. B 47, 1742 (1993). 
[23] W. P. Su, J. R. Schriefler and A. J. Heeger, Phys. Rev. 

Lett. 44, 1698 (1979); Phys. Rev. B 22, 2099 (1980). 
[24] H. Takayama, Y.R. Lin-Liu and K. Maki, Phys. Rev. B 

21, 2388 (1980). 



